图1 研究区地理位置图
下载摘 要：作者以黑龙江流域中俄跨境地区为研究区，以0、5、10和15 ℃为基础温度，基于大气再分析数据产品ECMWF Re-Analysis 5中的逐小时2 m气温数据，通过时空筛选与尺度转换、有效温度始末日期甄别，以及区间数据批量统计等步骤，得到黑龙江流域中俄跨境地区0.25°×0.25°有效积温数据集（2002–2020）。将本数据产品（有效积温）与基于气象实测数据及其关联产品的有效积温进行相关性分析表明，本产品与两者均呈正相关关系，并且年际间波动情况基本一致。本数据集内容包括：（1）研究区边界数据；（2）研究区2002–2020逐年0.25°×0.25°有效温度(≥0、≥5、≥10和≥15℃)始末时间数据；（3）研究区2002–2020逐年0.25°×0.25°有效积温(≥0、≥5、≥10和≥15℃)数据；（4）有效积温的验证数据。数据集存储为.shp、.tif和.xlsx格式，由305个文件组成，数据量为19.9 MB（压缩为1个文件，1.69 MB）。
关键词：黑龙江流域；中俄跨境地区；有效温度始末日期；有效积温；2002–2020
DOI: https://doi.org/10.3974/geodp.2022.03.18
CSTR: https://cstr.escience.org.cn/CSTR:20146.14.2022.03.18
数据可用性声明：本文关联实体数据集已在《全球变化数据仓储电子杂志（中英文）》出版，可获取：
https://doi.org/10.3974/geodb.2022.04.05.V1 或 https://cstr.escience.org.cn/CSTR:20146.11.2022.04.05.V1.
长期以来，农业生产的进程受到多种自然气候环境条件的影响，其中气温是作物发育的关键因素，其发生的任何变化都会影响地区气候资源的供热条件。气温与太阳辐射、光照周期和土壤温度等变量直接或间接相关，同时对作物光合作用过程起主导作用，所以气温也是影响作物生产力和质量的重要变量之一。除此之外，温度变化及其波动效应对农田生态系统研究也起到重要支撑作用，因其在影响农业资源供热变化的同时，还会导致火灾、霜冻、干旱等短时或持续性次生灾害发生，对整个生态系统的可持续性产生较大影响[1,2]。
逐日气温的累积效应是作物生长和生产潜力的一个重要指标[3]。大约300年前，de Réaumur提出了植物生长发育所需的蓄热统计理念，并由此产生一系列基于气温的评价指数来辅助对于各区域农业进行科学管理[4,5]。上述指数旨在将历史温度数据进行统计分析，以帮助人们了解各地区作物生长条件的差异，代表性方法包括平均生长季节温度、Huglin光热指数、活动累积温度和有效积温等方法。其中平均生长季节温度是将某一地区生长季内（北半球设为每年4–10月，南半球设为10月至下一年4月）的温度进行平均，进而形成对于当地农业气候适宜性的宏观认识[6]。Huglin光热指数是以10 ℃气温为基准评估葡萄种植栽培和酿酒所需的气候条件[7]。而活动累积温度是指日平均气温高于一定温度阈值（通常为10℃）的期间内高于零度气温的年累积总和，该指标较为直观地表征了某地的总体热量资源，因此其对于农业气候分析和区划等研究工作具有重要意义[8,9]。而有效积温法基于植物生长季节热量积累的理念，描述了高于生物学零度的日期内适宜其生长发育的日平均气温与其差值的年度累积情况，其中生物学零度代表了植物开始生长发育的最低要求[10]。在以上各种表示热量累积的评价方法中，有效积温法去除了地区热量资源中不适于作物生长的“无效”部分，进一步修正了两者间的线性关系，因此可用于监测作物生长条件的时空演变过程以及预测其未来生育时段，并可根据不同作物的生长温度阈值将相应的气温范围纳入区域乃至全球热势能的研究中，因而被广泛使用[11,12]。
黑龙江流域是横跨中国、俄罗斯、蒙古和朝鲜四国的一个跨境流域，其中位于黑龙江-乌苏里江沿线的中俄边境毗邻地区包含黑龙江省、阿穆尔州、犹太自治州以及滨海边疆区等行政单元（图1）。作为全流域的主要农业产区，该地自然资源储量较高、地理及气候条件总体良好，具有深厚的农业资源发展潜力。黑龙江省为温带大陆性季风气候，近年来平均气温为4 ℃左右，平均降水量范围介于530–760 mm之间，各类气象条件具有鲜明的季节性特征[13]。就其地理条件而言，松嫩与三江两大平原区位于全省东西两侧，内部土壤质量肥沃，适宜种植大豆、水稻、玉米、马铃薯等农作物。而受制于寒冷气温和永久冻土层等环境本底影响，俄罗斯远东区种植业主要集中于自然条件较好的南部和东部行政区，这些地区兼具温带大陆性和海洋性季风气候，大部分区域海拔低于500 m，其中泽雅-布列亚及兴凯湖平原为远东农耕重点地带，其雨热条件适宜种植谷类（大麦、燕麦、荞麦）、大豆、水稻等农作物，其中犹太自治州还大力发展果类种植[14]。从研究区内中俄跨境农业发展现状来看，两地近年来不断深化农业生产互补性合作，其中黑龙江省对俄农业合作种植面积超过950万亩，在远东从事农业开发的企业超过200家。双方科研单位的农业技术性交流以及农作物优势品种的进出口情况也呈上升态势，其中黑龙江省对阿穆尔州主产的大豆进口逐年大幅增长，有效缓解了我国对于该类作物的需求缺口[15]。而在当今气候变化的背景下，中俄两地的气温等农业气象资源条件不断变化，农业自然灾害的频率和严重程度也随之增加，对中俄跨境地区农业方面的贸易合作及农产品安全产生重大影响。故此有必要摸清该区域农作物生长发育的热量资源情况，从所处气候条件角度对于研究区农业生产活动、灾害防控及后备耕地潜力分析提供相关基础数据和理论参考。
图1 研究区地理位置图
下载《黑龙江流域中俄跨境地区0.25°×0.25°有效积温数据集（2002-2020）》[16]的名称、作者、地理区域、数据年代、时间分辨率、空间分辨率、数据集组成、数据出版与共享服务平台、数据共享政策等信息见表1。
| 条目 | 描述 |
|---|---|
| 数据集名称 | 黑龙江流域中俄跨境地区0.25°×0.25°有效积温数据集（2002-2020） |
| 数据集短名 | EffecAccTemp_CR_2002-2020 |
| 作者信息 | 周业智0000-0002-6826-151X，中国科学院地理科学与资源研究所，zhouyz@lreis.ac.cn 王卷乐0000-0002-5641-0813，中国科学院地理科学与资源研究所，wangjl@igsnrr.ac.cn 李凯0000-0001-6234-6806，中国科学院地理科学与资源研究所，lk@lreis.ac.cn Grigorieva, E. A. 0000-0002-7811-7853，俄罗斯科学院区域问题复杂分析研究所远东地区分所，eagrigor@yandex.ru |
| 地理区域 | 黑龙江流域内部中俄跨边境行政区域（黑龙江全省、犹太自治州、阿穆尔州和滨海边疆区部分地区） |
| 数据年代 | 2002–2020 |
| 时间分辨率 | 1年 |
| 空间分辨率 | 0.25°×0.25° |
| 数据格式 | .shp、.tif、.xlsx |
| 数据量 | 19.9 MB（压缩后） |
| 数据集组成 | 本数据集内容包括：（1）研究区边界数据；（2）研究区2002–2020逐年0.25°×0.25°有效温度(≥0℃、≥5℃、≥10℃和≥15℃)始末时间数据；（3）研究区2002–2020逐年0.25°×0.25°有效积温(≥0℃、≥5℃、≥10℃和≥15℃)数据；（4）有效积温的验证数据 |
| 基金项目 | 中国科学院（XDA2003020302）；中国工程院（CKCEST-2022-1-41） |
| 数据计算环境 | Python，ArcGIS |
| 出版与共享服务平台 | 全球变化科学研究数据出版系统 http://www.geodoi.ac.cn |
| 地址 | 北京市朝阳区大屯路甲11号100101，中国科学院地理科学与资源研究所 |
| 数据共享政策 | （1）“数据”以最便利的方式通过互联网系统免费向全社会开放，用户免费浏览、免费下载；（2）最终用户使用“数据”需要按照引用格式在参考文献或适当的位置标注数据来源；（3）增值服务用户或以任何形式散发和传播（包括通过计算机服务器）“数据”的用户需要与《全球变化数据学报》（中英文）编辑部签署书面协议，获得许可；（4）摘取“数据”中的部分记录创作新数据的作者需要遵循10%引用原则，即从本数据集中摘取的数据记录少于新数据集总记录量的10%，同时需要对摘取的数据记录标注数据来源[17] |
| 数据和论文检索系统 | DOI，CSTR，Crossref，DCI，CSCD，CNKI，SciEngine，WDS/ISC，GEOSS |
制作本数据集的原始气象数据产品来源于欧洲中期天气预报中心[1]发布的全球逐小时2 m空气温度数据，数据时间为2002年至2020年。结果验证数据来源于基于CN05.1格点数据资料计算的2002-2018年中国黑龙江省对应有效积温数据以及Elena Grigorieva基于气象站点观测的逐日气温数据制作的2002-2011年俄罗斯远东南部地区有效积温数据。
有效积温（effective accumulated temperature，EAT）是与植物生长息息相关的评价指数，其通过统计高于生物学零度的日平均气温（有效温度）累积值来反映供目标植物适宜生长的热量资源储量。该指数的构建过程包含了植物在生长季节所经历的时间和温度等信息，故此该指数的变化将对植物物候期产生影响，进而导致各区域农业生产活动及其结果的变化。具体组合计算公式如式1、式2所示，而现有研究表明，利用小时级记录的温度数据进行逐日平均温度计算将会更为准确的代表某地某天的气温状况，从而对于该要素的计算结果进行更精确的表征[18,19]。故此本研究中关于日平均气温的计算方法如式3所示：
(1)
(2)
(3)
式中，Ti为按不同统计方法计算的第i天的日平均温度（℃），Tbase为目标作物适宜生长的基础温度（℃）。i =1, 2, …, m为“逐日平均气温≥生长基础温度”的时段范围，Tmax和Tmin为当天的最高和最低温度（℃），Thour为一天内逐小时记录的气温数据。
图2为本研究的整体工作流程，由数据获取与预处理、作物生长界限温度选定、有效温度始末时间确定与统计以及成果质量检验等四部分组成，具体工作流程如下：
图2 数据研发技术流程图
下载在制作本数据集的原始数据处理过程中，本研究利用python语言和ArcGIS平台基于研究区矢量范围完成对研究时段内ERA-5数据的批量裁剪，而后将一天内的逐小时温度数据进行平均得到2002–2020年该区域逐日平均气温。同时本研究对所采用的验证数据按类型进行栅格数据裁剪及面板数据筛选等预处理操作。
植物会随着长时间的推移显现出其对于各量级温度的适应能力，当其周围温度高于一定数值时会开始发育，反之在其以下则鲜有生长。经过相关研究者对各类型作物的观测实验与数据整理，其中主要经济作物的生长界限温度如表2所示：
| 作物类型 | 生长界限温度Tbase（℃） |
|---|---|
| 卷心菜 | 0 |
| 豌豆 | 4.5 |
| 小麦 | 5 |
| 马铃薯 | 7 |
| 玉米 | 10 |
| 大豆 | 10 |
| 水稻 | 15 |
结合研究区实际种植的作物类型，本文最终选定Tbase=0、5、10和15 ℃来表征研究区多种有效积温的起始（结束）的“基线”温度。
在某一年有效积温始末时期的确定流程中，本文首先根据研究区的历史物候和生育期统计信息将有效温度开始与结束日期出现的时间段分别设定为每年1–7月和8–12月。然后选择5日滑动平均值（该年内每天取当日、前2日和后2日的5日序列计算平均值）大于设定的界限温度的最长一段时期作为有效温度起始（终止）范围，在起始范围的首个序列中，挑选第一个Ti大于界限温度的日期作为起始日；在终止范围的最后一个序列中，选取最后一个Ti大于界限温度的日期作为终止日。
在有效积温的计算过程中，除考虑作物的生长基础温度Tbase外，其在生育期内的最高上限温度TUT也会对作物生长情况产生影响。先前研究显示极端高温造成的热浪天气会对各类作物的生长发育造成一定的减缓与抑制[21,22]。综合考量研究区主要作物各发育阶段的试验资料及相关文献[11,20,23]，本研究选定Ti=30 ℃作为研究区有效积温统计的上限温度。
在Ti大于Tbase的时间范围内，将小于30 ℃的Ti与四种Tbase的差值进行累加获得研究区逐年有效积温（图3）。后将积温统计结果与和基于实测数据制作的验证结果进行相关性分析，以评估本数据产品对该区域热量资源表征的质量和效果。
图3 本研究有效积温计算原理示意图
下载本数据集的边界矢量数据存储为Shapefile文件，数据名称为studyarea.shp。有效温度始末时期及有效积温数据存储为单波段Geotiff文件，其中四种有效温度开始时间的数据名称格式为atTbasefirstday_EAT_0.25°_YYYY.tif；结束时间的数据名称格式为atTbaselastday_EAT_0.25°_YYYY.tif；有效积温数据名称格式为atTbase_EAT_0.25°_YYYY.tif。其中Tbase表示本研究设定的四种生长基础温度，EAT表示数据内容为有效积温结果；0.25°表示数据集的空间分辨率为0.25°×0.25°；YYYY表示数据文件所处时间为YYYY年，数据集其他信息如表3所示。验证数据中，黑龙江省有效积温数据存储为单波段Geotiff文件，名称格式与研究数据集中的对应结果相同，而基于气象站点统计的俄罗斯远东南部地区有效积温数据集存储为.xlsx文件。
| 编号 | 属性 | 研究数据数值 |
|---|---|---|
| 1 | 像素类型 | 浮点型 |
| 2 | 行数 | 67 |
| 3 | 列数 | 50 |
| 4 | 像素值 | 0-4,000 |
| 5 | 像素大小 | 0.25°×0.25° |
| 6 | 坐标系统 | WGS84 |
通过先前的数据处理流程得到研究区四类有效积温年度空间分布及趋势变化情况，其中≥0 ℃的情形下的对应结果如图4所示。
图4 研究区2002–2020年≥0℃的有效积温空间分布（4a-4s）及其历年变化趋势图（4t）
下载由上图可知，研究区内中俄两地≥0 ℃的有效积温空间格局总体呈现由北向南增加的态势，其中大于3,000 ℃的积温高值区主要分布于黑龙江省东西两侧、俄罗斯远东阿穆尔州南部以及滨海边疆区西部等地，与研究区耕地分布情况总体吻合。低值区位于阿穆尔州西北部高海拔地区，该地常年积温位于1,000 ℃以下。而就19年间的积温变化趋势而言，研究区大部分区域积温变化呈现增长态势，其中幅度较大（>20 ℃/4a）的区域出现在阿穆尔州南部、犹太自治州北部以及黑龙江省与滨海边疆区毗邻地区。而黑龙江北部与阿穆尔州交界以及该省西部地区年际积温较为稳定，范围内部分地区呈现出轻微下降的趋势（5–10 ℃/4a）。
从下图5a、5b中可以看出，研究区≥0 ℃有效温度的最早开始和结束的时间及地点分别为3月和9月中旬，对应的最晚时间分别为5月和11月上旬。其中阿穆尔州东北部区域积温时间最短，而黑龙江省南部区域积温时间相对较长，持续时间最长的地区位于该省西南及东部与俄罗斯毗连地带。就研究区近20年前后两时段的对比情况可以发现，在黑龙江省东西两侧、俄罗斯犹太自治州和滨海边疆区西部的农业种植区近10年内的平均积温持续时间较前一时段延长（开始时间提前，结束时间推迟），而俄罗斯阿穆尔州南部对应区域的积温持续时间较前期则有所缩短。
图5 ≥0℃和≥10℃的研究区2011–2020年平均有效积温开始（5a、5c）与结束时间（5b、5d）空间分布（黑色和红色散点分别代表对应位置较2002–2010年平均结果的提前与推迟）
下载研究区近10年间>10℃有效温度最早开始于4月中旬，最晚结束于10月上旬。黑龙江省大部分地势平坦地区的有效温度起始时间均早于5月份，而由西北向东南走向的高海拔区域的对应时间普遍延后10天左右。俄罗斯远东区域内的对应日期基本始于五月之后，阿穆尔州南部和滨海边疆区以西部分农业种植区积温的起始时间则略有提前。在积温结束时间的分布情况中，黑龙江南部、犹太自治州和滨海边疆区大部普遍晚于九月末，而阿穆尔州除其与黑龙江省毗邻区域外的积温时间结束于九月上旬左右。对比本研究前后两时段的始末时间可以发现，全区绝大部分地区近10年间起始时间较前期提前，而在结束时间的对比结果中，北部的俄罗斯两州出现推迟情况的样点占比超过90%，而黑龙江省与滨海边疆区两地积温末期的提前与延后比例基本持平，其中出现提前的地区主要位于黑龙江省以西和东部中俄跨境区域。
根据本研究中≥10 ℃的有效积温热量标准，可将研究区逐年动态的划分为四个积温区间（对应at 10℃_EAT的范围为>1,400 ℃，1,400–1,100 ℃，1,100–800 ℃及<800 ℃），并对其多年面积占比变化情况进行统计与分析（图6）：
图6 基于≥10℃有效积温标准的各区划年际所占面积比例
下载由上图可知，研究区四类积温区划面积占比大致为7.5%，27.5%，35%和30%。其中一、三积温区变化趋势较为稳定，而第二（四）积温区年际变化呈现略增长（下降）态势。在各区划年际间波动较为明显的时间点中，2010年积温高值区面积占比较高，其在第一、二积温区的面积分别超出15%和35%。而与之相反的情况出现在与之相邻的2009年，当年有效积温处在1,100 ℃以下的面积比例约占整个研究区的7成左右。
在数据产品的验证过程中，研究基于国家气候中心制作的CN05.1数据产品中所包含的逐日平均气温按上述处理流程对黑龙江省的2002–2018年有效积温进行计算并与本研究对应区域的结果进行对照分析。CN05.1数据是基于1961–2018年间中国大陆范围内2,400余个国家级气象台站实测数据通过薄盘样条函数外推插值法制作的产品。先前研究得出的积温与中国生态系统研究网络中积温站统计的对应结果具有较好的相关性[24]。俄罗斯远东地区的验证数据采用Elena Grigorieva基于气象站点观测的逐日气温数据制作的2002–2011年俄罗斯远东南部地区四种对应有效积温数据[11]。在验证前的数据预处理过程中，本研究将计算得到研究区及黑龙江省的栅格数据产品基于各气象站点位置通过双线性插值法进行数值提取，而后连同俄罗斯远东部分的计算结果进行相关性分析，对应结果如图7、图8所示。
图7 ≥0℃有效积温的评估结果
下载
图8 ≥10℃有效积温的评估结果
下载由以上两类有效积温的对比验证结果可知，在2002–2020年间，除个别站点的数据有较大差异外，本研究计算得到的≥0 ℃和≥10 ℃的有效积温与基于CN05.1和实测温度数据得出的对应结果的多年平均值的体量基本一致（图7a、8a）。而数据间的平均误差（Mean Error，ME）结果（图7b，8b）显示，本研究在黑龙江省内大部分站点处的计算结果较对照数据存在高估，而在俄罗斯远东地区的情况则相反。在两类积温数据中，除黑龙江省第7、8、19和23号站点以及俄罗斯远东第3和6号站点外，剩余37座站点间的ME值均落入-200–200 ℃的范围。从两类积温数据组内的Pearson相关系数结果（图7c、8c）来看，大部分站点的r值在0.45以上，基本反映出近20年来研究区两类积温年际间的波动变化状态。在各组对应的均方根误差（Root Mean Squared Error，RMSE）统计结果（图7d、8d）中，两组数据的最大偏差范围未超过500 ℃，且绝大部分站点的两类RMSE值处于各自最大值的50%以内。以上条件说明本研究与以上两种基于实测数据计算的积温结果拟合程度良好，可作为反映年尺度下该区域不同作物所需热量资源的基础数据。
本研究首先针对ERA5再分析产品中的逐小时2 m温度数据进行加工转换，将得到的逐日平均温度与设定的作物生长发育基础温度（0、5、10和15 ℃）按照对应计算原理对黑龙江流域中俄跨境地区2002–2020年间有效温度产生的始末时间进行确定并统计各生育期的年度有效积温，进而从地物类型和地貌条件等角度系统分析了两者的空间分布特征，并对各级有效积温区划面积在年际间的动态演变规律进行了宏观刻画。结果表明有效积温高值主要位于研究区南部平原地区，与研究区农耕地分布情况基本一致。而近20年内研究区大部分地区积温变化呈现稳定增长趋势，其中西北部地区在前三种情形下的对应情况更为明显，未来可就这一条件结合其他相关评价因素对于此类非农耕地域的后备潜力进行分析挖掘，进而加强在气候变化和国际局势等不确定条件下对于中俄跨边境农业发展和粮食安全的保障能力。
作者分工：王卷乐负责数据集开发总体设计工作；周业智、李凯、Grigorieva E.A负责数据采集、处理及验证工作；王卷乐、周业智、李凯共同完成数据论文撰写。
致谢：感谢中国科学院特别交流计划项目对本研究的支持。
利益冲突声明：本研究不存在研究者以及与公开研究成果有关的利益冲突。
[1] European Centre for Medium-Range Weather Forecasts. https://www.ecmwf.int/en/forecasts/dataset/ecmwf-reanalysis-v5.